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ABSTRACT 

Turbulent heating may play an important role in galaxy-cluster plasmas, but 
if turbulent heating is to balance radiative cooling in a quasi-steady state, some 
mechanism must regulate the turbulent velocity so that it is neither too large 
nor too small. This paper explores one possible regulating mechanism associ- 
ated with an active galactic nucleus at cluster center. A steady-state model for 
the intracluster medium is presented in which radiative cooling is balanced by 
a combination of turbulent heating and thermal conduction. The turbulence is 
generated by convection driven by the buoyancy of cosmic rays produced by a 
central radio source. The cosmic-ray luminosity is powered by the accretion of in- 
tracluster plasma onto a central black hole. The model makes the rather extreme 
assumption that the cosmic rays and thermal plasma are completely mixed. Al- 
though the intracluster medium is convectively unstable near cluster center in 
the model solutions, the specific entropy of the thermal plasma still increases 
outwards because of the cosmic-ray modification to the stability criterion. The 
model provides a self-consistent calculation of the turbulent velocity as a function 
of position, but fails to reproduce the steep central density profiles observed in 
clusters. The principal difficulty is that in order for the fully mixed intracluster 
medium to become convectively unstable, the cosmic-ray pressure must become 
comparable to or greater than the thermal pressure within the convective re- 
gion. The large cosmic-ray pressure gradient then provides much of the support 
against gravity, reducing the thermal pressure gradient near cluster center and 
decreasing the central plasma density gradient. A more realistic AGN- feedback 
model of intracluster turbulence in which relativistic and thermal plasmas are 
only partially mixed may have greater success. 
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1. Introduction 

Recent x-ray observations show that some source of heating inhibits cooling of intra- 
cluster plasma to temperatures below roughly one-third of the average cluster temperature 
(Peterson et al 2001, Peterson et al 2003, Bohringer et al 2001, Tamura et al 2001, Molendi 
& Pizzolato 2001). One of the relevant heating mechanisms is conduction of heat from a 
cluster's hot outer region towards cluster center. Fabian et al (2002), Voigt et al (2002), and 
Zakamska & Narayan (2003) showed that thermal conduction can balance radiative cooling 
in many clusters if the thermal conductivity k,t is in the range 0.1 — 1.0/ts, where k,s is 
the Spitzer thermal conductivity. However, models in which cooling is balanced entirely by 
conduction face some difficulties. Different clusters need different values of k,t to fit the 
observations, and it is not clear why kt would vary in the required way. Relatively cool 
cluster plasmas can require values of kt in excess of kq, for conduction to balance cooling, 
values which are difficult if not impossible to justify. Theoretical studies of heat transport 
in tangled intracluster magnetic fields find values of k t towards the lower end of the range 
of required values (Narayan & Medvedev 2001, Chandran & Maron 2004, Bessho, Maron, & 
Chandran 2004, and references therein). In addition, pure conduction models are thermally 
unstable, although the implications for clusters are debated since the instability growth time 
is long (~ 2 - 5 Gyr for k t ~ 0.2 - 0Ak s ) (Kim & Narayan 2003, Soker 2003). 

Other heating mechanisms include dissipation of turbulent motions (Loewenstein & 
Fabian 1990, Churazov et al 2004, Dennis & Chandran 2004) and turbulent mixing (Cho 
et al 2003, Voigt k Fabian 2004, Kim & Narayan 2003, Dennis & Chandran 2004, Narayan 
& Kim 2004). Observational evidence of plasma motion in the Perseus cluster at roughly 
half the sound speed c s (Churazov et al 2004) provides support for models invoking turbulent 
heating, but also poses a theoretical challenge for the following reason. As the rms turbulent 
velocity u varies from 0.1c s to c s , and as the dominant velocity length scale at radius r varies 
from O.lr to r, the turbulent heating rate in a typical cluster varies from a value much smaller 
than the radiative cooling rate it! to a value much larger than R (Dennis & Chandran 2004). 
If turbulent heating in fact balances radiative cooling, an explanation is needed for how the 
turbulence amplitude is fine-tuned to the required value. 

A number of authors have considered plasma heating by an active galactic nucleus 
(AGN) at a cluster's center. Several mechanisms have been considered for transferring 
AGN power to the ambient plasma, including buoyantly rising bubbles of cosmic rays or 
centrally heated gas (Tabor & Binney 1993, Churazov et al 2000, Churazov et al 2001, 
Begelman 2001a, 2001b, Reynolds 2001, Bruggen et al 2002, Begelman 2002, Ruszkowski 
& Begelman 2002, Fabian et al 2003, Mathews et al 2003, Reynolds et al 2004), wave- 
mediated plasma heating by cosmic-ray electrons (Rosner & Tucker 1989) or cosmic-ray 
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protons (Bohringer & Morfill 1988, Loewenstein, Zweibel, & Begelman 1991), turbulence 
(Loewenstein & Fabian 1990, Churazov et al 2004), Compton heating (Binney & Tabor 
1995, Ciotti & Ostriker 1997, 2001), and jet mechanical luminosity and shocks (Binney & 
Tabor 1995). An attractive feature of AGN- heating models is that the heating rate increases 
with the mass accretion rate of intracluster plasma, which allows the models to produce 
globally stable equilibria. If the mass accretion rate rises above the equilibrium value, AGN 
heating increases, and radiative cooling of the ICM is less able to cause the net cooling that 
drives mass accretion onto the central supermassive black hole. If the mass accretion rate 
falls below the equilibrium value, the net cooling rate rises, restoring the mass accretion 
rate to its equilibrium level. If the thermal conductivity is not dramatically reduced relative 
to the Spitzer value, small-scale thermal instabilities are also suppressed (Rosner & Tucker 
1989, Ruszkowski & Begelman 2002, Zakamska & Narayan 2003). A powerful motivation for 
AGN-heating models is the observation that almost all clusters with strongly cooling cores 
possess active central radio sources (Eilek 2003). 

The present paper builds upon previous studies of turbulent heating and AGN feedback 
in clusters, and explores one possible explanation for how the turbulent velocity could achieve 
the value needed to balance radiative cooling. A steady-state model of the ICM is presented 
in which radiative cooling causes inflow of intracluster plasma and accretion onto a central 
supermassive black hole, leading to cosmic-ray production by a central radio source. The 
cosmic-ray luminosity L cr is taken to be proportional to the mass accretion rate M, and 
it is assumed that the cosmic-rays mix completely into the thermal intracluster plasma. 1 
Since the cosmic rays provide pressure without noticeably increasing the density of the ICM, 
they increase the buoyancy of the plasma, eventually leading to convection. Convection 
heats the plasma in three ways, through the viscous dissipation of turbulent motions, by 
mixing hot plasma from the outer regions of a cluster in towards cluster center, and by 
providing a vehicle for cosmic-ray pressure to do work on the thermal plasma. Convection 
is treated with a two-fluid (thermal plasma and cosmic ray) mixing length theory developed 
in section 2. The turbulent velocity and turbulent heating rate increase with increasing M 
and L cr . In equilibrium, M and L CT attain those values for which turbulent heating and 
thermal conduction balance radiative cooling. This equilibrium is expected to be stable for 
the same reasons as other AGN- feedback/thermal-conduction models. 

Once the outer temperature and density are specified, the model can be used to calculate 
the density, temperature, and cosmic-ray pressure throughout the cluster. The model density 
profiles, however, are not as steep as observed profiles within the central ~ 50 kpc, where 



See Ensslin (2003) for a detailed discussion of the escape of cosmic rays from radio-galaxy cocoons. 
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the ICM is convective. The reason for this discrepancy is that the cosmic-ray pressure has 
to be comparable to or greater than the thermal pressure within the convective region in 
order for cosmic-ray buoyancy to make the ICM convectively unstable. The large cosmic-ray 
pressure gradient provides much of the support of the ICM against gravity, reducing the 
thermal pressure gradient, and decreasing the central plasma density. A similar discrepancy 
was found in an MHD- wave-mediated cosmic-ray heating model (Loewenstein et al 1991), 
which also involved significant nonthermal pressure. The present model also provides a self- 
consistent calculation of the profiles of the turbulent velocity and turbulent heating rates in 
the ICM. The resulting velocities are subsonic, implying that convective regions are fairly 
close to marginal convective stability. However, because the cosmic-ray pressure modifies 
the stability criterion, the specific entropy of thermal plasma still increases outwards. 

The discrepancy between the model and observations may be linked to the assumption 
that cosmic rays are completely mixed into the thermal plasma. This assumption, which 
is made to simplify the analysis, is inconsistent with the x-ray cavities (depressions in x- 
ray emission often associated with enhanced synchrotron emission) observed in roughly one- 
fourth of the clusters in the Chandra archive (Birzan et al 2004). When mixing is incomplete, 
cosmic-ray pressure is less able to support the thermal plasma, the thermal pressure gradient 
has to be larger, and therefore the central plasma density has to be larger. A model that 
accounts for incomplete mixing may have greater success. 

The remainder of the paper is organized as follows. In section 2, the model is de- 
scribed in detail and example solutions are compared to observations of Abell 478. Section 3 
summarizes the results of the paper. 



In the cooling flow model of intracluster plasmas (Fabian 1994), entropy-generating 
heat sources are neglected, and radiative losses are balanced by inward advection of plasma 
internal energy as well as gravitational and pdV work. The resulting average radial velocity 
corresponds to mass accretion rates above 1OOOM yr _1 in some clusters. In contrast, in 
the present model, radiative cooling is balanced by heating from turbulence and thermal 
conduction, leading to much smaller radial velocities and mass accretion rates. The average 
fluid velocity is neglected in the fluid equations, and the mass accretion rate is set equal to 
Bondi (1952) rate corresponding to the equilibrium plasma parameters at cluster center, 



2. 



The two-fluid convection model 



M = 



7rG 2 M b V(0) 

C s ,ad(0) 3 



(1) 
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where M hh is the mass of a supermassive black hole at cluster center, p(0) is the central 
density, and c S;a( j(0) is the adiabatic sound speed at cluster center. 2 

The inferred mass accretion rate is used to determine the cosmic-ray luminosity through 
the equation 

L CT = V Mc 2 , (2) 

where 77 is a dimensionless efficiency factor. The spatial distribution of cosmic-ray injection 
into the ICM is a major uncertainty in the model. Some clues are provided by radio ob- 
servations, which show that cluster-center radio sources (CCRS) differ morphologically from 
radio sources in other environments. As discussed by Eilek (2003), roughly half of the CCRS 
in a sample of 250 sources studied by Owen & Ledlow (1997) are "amorphous," or quasi- 
isotropic, presumably due to jet disruption by the comparatively high-pressure cluster-core 
plasma. With the exception of Hydra A, the CCRS in the Owen-Ledlow (1997) study are 
smaller than non-cluster-center sources, with most extending less than 50 kpc from cluster 
center (Eilek 2003). Complicating matters is the possibility that cosmic-ray bubbles rise 
buoyantly away from their acceleration site before mixing into the thermal plasma, effec- 
tively distributing cosmic-ray injection over a larger volume. In this paper, the cosmic-ray 
energy introduced into the intracluster medium per unit volume per unit time is taken to be 

S(r) = S e- r2 ' r *, (3) 

where the constant r s is a free parameter and the constant So is determined from the equation 
L CT = 4n J °° drr 2 S(r) and equation (2). 

The magnetic pressure is taken to be much less than the thermal pressure. Magnetic 
terms are then neglected in the fluid equations, and the effects of the magnetic field on con- 
vection are ignored. 3 Although the magnetic field does not affect the bulk fluid velocity, for 
any plausible field strength the field plays an important role in heat transport by constrain- 
ing charged particles to move primarily along field lines. Since the magnetic field in clusters 
is disordered and probably turbulent (Kronberg 1994, Taylor et al 2001, 2002), chaotic field- 
line trajectories inhibit large radial particle excursions, suppressing radial heat transport to 
some degree. In this paper, the thermal conductivity is set equal to the Spitzer (1962) value 
for a non- magnetized plasma, Kg, multiplied by a suppression factor 9, 

k t = 6k s , (4) 



2 See Quataert & Narayan (2000), Nulsen (2003), and Bohringer et al (2003) for further discussions of 
Bondi accretion in clusters. 

3 See Pen et al (2003) for a very different picture of the interplay between magnetic fields and convection 
in quasi-spherical accretion, in which convective motions are resisted by magnetic tension. 
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where 



n e /c# \5 keV/ \ n e J \lnA c 

n e is the electron density, T is the temperature, k B is the Boltzmann constant, and In A c is the 
Coulomb logarithm. Recent theoretical studies suggest that 9 ~ 0.2 (Narayan & Medvedev 
2001, Chandran & Maron 2004, Maron, Chandran, & Blackman 2004). The heating rate per 
unit volume from thermal conduction is 



H tc = V-(ktVT). (6) 

Following Tozzi & Norman (2001) and Ruszkowski & Begelman (2002), the radiative cooling 
per unit volume per unit time for free-free and line emission is taken to be 

• 10~ 22 ergs cm 3 s"\ (7) 

where n,i is the ion density and the numerical constants correspond to 30% solar metallicity. 
Radiative cooling of cosmic rays is ignored, which is reasonable if protons make the dominant 
contribution to the cosmic-ray pressure. Cosmic rays are assumed to diffuse relative to the 
thermal plasma, but the value of the cosmic-ray diffusion coefficient D CI in clusters is not 
known. In the model, 

Dcr = ^D$ + vy, (8) 

where D and Vd are constants. 

The thermal plasma and cosmic rays are treated as co-extensive fluids that interact 
only in limited ways. Wave pitch-angle scattering causes the two fluids to move with the 
same bulk fluid velocity v, although the finite efficiency of scattering allows for cosmic-ray 
diffusion as already discussed. Collisional and microphysical collisionless energy exchange 
between the two fluids is neglected, but the cosmic rays can heat the plasma, and vice versa, 
through pdV work. The thermal plasma is treated as a fluid of adiabatic index 7, 

V = (7 - l)pe, (9) 

where p is the plasma pressure, p is the plasma density, and e is the plasma internal energy 
per unit mass. The cosmic rays are treated as a fluid of adiabatic index 7 cr , 

Per = (lev - l)Pcre C r, (10) 

where p cr is the cosmic-ray pressure, p CT is the cosmic-ray density, and e cr is the cosmic- 
ray internal energy per unit mass. Both 7 and 7 cr are treated as constants (although the 
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0.0086 
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-1.7 



+ 0.058 



k B T 
1 keV 



0.5 



+ 0.063 
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adiabatic index of a fluid should really vary with temperature from 5/3 in the non-relativistic 
limit to 4/3 in the ultra- relativistic limit). 



The first law of thermodynamics can be written 

p dp 



de 
dt 



H + 



p 2 dV 



(11) 



where d/dt = (d/dt + v • V), and H is the net heating per unit mass per unit time. 
Equations (9) and (11) together give 



1 



7 



dp 7p dp 
dt p dt 



pH. 



The equation of continuity is 



dp 
dt 



+ V • (pv) = S M , 



(12) 



(13) 



where Sm is the rate at which mass is introduced (by, e.g., stellar winds) per unit time per 
unit volume. Equations (12) and (13) give 



7 



dp 
dt 



+ V • (vp) 



-pV • v + pH + 



-fpS, 



M 



Similarly, 



dp a 

dt 



+ V • (vp cr ) 



-PcrV • V + p cr H cr + 



(t-i)p' 

IcvPcvSm X 



(14) 



(15) 



7 cr - 1 [ dt (7 cr - l)p cr ' 

where H CT and Sm,ct are, respectively, the cosmic-ray heating and mass-injection rates. For 
the thermal plasma, 

pH = H diss + H tc - R, (16) 

where i/diss is the rate of heating from viscous dissipation of turbulent motions, quantified 
below in equation (50). For simplicity, it is assumed that 



S 



M 



0. 



The right-hand side of equation (15) is simplified by assuming that 

7crPcr>S'M,( 



(7cr - l)Pc 



= V-(£> cr Vp cr )+S(r), 



(17) 



(18) 



with D cr given by equation (8) and S(r) given by equation (3). The first term on the right- 
hand side of equation (18) is a simplified representation of diffusive energy transport, and 
the second term models the energy injection from the central radio source. 
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Convection and convective stability are treated using the following simple two-fluid 
mixing length theory. Consider a parcel of plasma and cosmic rays initially at a distance r = 
ro from cluster center. The parcel is displaced radially by an amount I (positive or negative) 
to r = ri = r + Z, where |/| is the mixing length. For simplicity, thermal conduction and 
cosmic-ray diffusion into and out of the parcel are neglected at this stage, so that the parcel 
expands adiabatically. The parcel is then mixed with its surroundings. The values of p, p, 
p cr , and p cr in the parcel are initially the same as the average values at r = r , denoted p , 
Po, Pct,o and p CT fi- After the parcel is displaced radially outwards, the new fluid quantities 
within the parcel are denoted p', p', p' cr , and p' CI . The average fluid quantities at r = r\ are 
denoted pi, pi, p cr ,i, an d p C r,i- The difference between the thermal-plasma density in the 
displaced parcel and its surroundings at r — T\ is denoted 

Ap = p'- Pl . (19) 

Similarly, Ap — p f — pi, etc. Since the parcel expands adiabatically, 

P' = Po(£f, (20) 

and 

Per = Pcr,0 ( • (21) 
VPcr,0/ 

The volume occupied by the thermal plasma and the volume occupied by the cosmic rays 
expand by the same amount, which implies that 



1 + 5. (22) 

PO Pcrfl 

It is assumed that the turbulent velocities are subsonic, so that the total pressure 

Ptot = P + Per (23) 

inside the parcel remains approximately the same as the total pressure outside the parcel. It is 
now assumed that |/| <r, which implies that 5 1. To lowest order in \l\/r, equations (20) 
through (23) give 

5 = —1 (24) 

7P + 7crPcr dr 

l P d Ptot jdp 

IP + lcrPcr dr dr' 

Ap _ hP dp tot ^dp 

7P + 7crPcr dr dr' 
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and 

Ap cr = -Ap. (27) 

Since only the lowest-order terms have been kept, it is not necessary to specify in equa- 
tions (24) through (26) whether p, p, etc are evaluated at r or ri, and so the subscripts on 
p, p, etc have been dropped. 

Since the displacement I is a signed quantity, the criterion for convective stability is 

^ > 0. (28) 
If p >> p CT , equation (28) is equivalent to the Schwarzchild criterion 

V(^>0. (29) 



dr 

More generally, the ICM can be convectively unstable even when the specific entropy of 
the thermal plasma increases outwards, provided that dp CT /dr is negative and \dp CT /dr\ is 
sufficiently large. 

The rms turbulent velocity u is assumed to satisfy 

2 ( *e£ if r'A„<o 

pu ~ < 8 (30) 
[ if l~ 1 Ap>0, 

where g is the gravitational acceleration. That is, the bulk-flow kinetic energy of a moving 
parcel in a convective region is approximately the buoyancy force on the fully displaced 
parcel times the mixing length times the standard numerical coefficient used in mixing- 
length theory (Cox & Giugli 1968). Equation (30) is then modified slightly so that du/dr 
varies continuously to zero as Ap increases through zero: 

„2< < + f if r'Ap<0 



u e 



where 



~ aq if l^Ap > 0, 



* = ^, (32) 
P 

° = -A, ( 33 ) 



8u 



u = u' r, (34) 
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and u' Q is a constant chosen so that uq remains much smaller than the sound speed throughout 
the cluster. The value of a in equation (33) is chosen so that du 2 /dq is continuous at q — 0. 

The approximate (signed) value of |V • v|, denoted (div v), is given by the fractional 
change in the volume of a parcel as it rises a distance I, [(p') 1 — Po 1 ]/(Po 1 )> divided by the 
time for the parcel to rise, \l/u\: 

(A - \ sgn(/)w dptot . . 

(div v) = — , (35) 

1P + 1ctPct dr 

where sgn(Z) = 

From equation (14), the steady-state thermal-plasma energy equation is 

• (vp) = -pV • v + H diss + H tc - R. (36) 
The steady-state cosmic-ray energy equation, from equations (15) and (18) is 

-V • (vp cr ) = - Pcr V • v + V • (D cr V Pcr ) + S. (37) 



7cr - 1 

To obtain averaged equations, each fluid quantity is written as an average value plus a 
turbulent fluctuation: 

v = (v) + Sv, (38) 

p = (p) + Sp, (39) 

etc. As mentioned above, (v) is set equal to zero. It is assumed that averaged quantities 
depend only on the radial coordinate r. The average (5v Sp), which is 7 — 1 times the 
thermal-plasma internal-energy flux, is estimated to be c m i x [w sgn(Z)f] Ap, where c m i x is a 
constant of order unity, giving 

(<Jv Sp) = rDQ, (40) 

where 

D = c mix u\l\ (41) 

is the eddy diffusivity, and 

= Ap = 7P dp tot _ dp 
I >yp + 7 cr p cr dr dr' 

The average (Sv Sp CT ) is estimated to be c m i x [u sgn(/)f] Aj9 cr , which yields 

(SvSp CI ) = -(5vSp). (43) 
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The average (SpV ■ 5v) is taken to be c mix Ap (div v), or 

(6pS7-6v) = DQ dp , tot . (44) 

7P + 7crPcr dr 

Similarly, (<5p cr V • 5v) is taken to be c mix Ap cr (div v), giving 

(5 Pa V -5v) = - (5 P V ■ 5v) . (45) 

Equations (43) and (45) reflect the assumption that the the total-pressure fluctuation in a 
displaced fluid parcel vanishes. The mixing length is set equal to 

\l\ = ar, (46) 

where a is a constant. Although |/| <C r was previously assumed, it is now assumed that a 
is of order unity, an inconsistency that also characterizes standard mixing-length theory. 

Upon averaging equations (36) and (37), discarding terms proportional to the average 
velocity, and dropping the angle brackets (...) around averaged quantities, one obtains 

Tr^-rX rD Q) = ; — + H d]ss + H tc - R, (47) 

(7 - l)r 2 dr jp + j CT p cr dr 

and 

-^i-l^^^^^AU^ + S. (48) 
(7 cr — l)r z dr IP + 7crPcr dr r z dr \ dr J 

In the averages of H tc and R the density and temperature are simply replaced by their average 
values. The average and turbulent velocities are neglected in the momentum equation, giving 

dPtot d$ , 4g ^ 

dr dr ' 

where $ is the gravitational potential, which is assumed to be dominated by a fixed dark 
matter distribution. The average value of H diss is set equal to 

#diss = j , (50) 

where Cdi SS is a dimensionless constant of order unity. Equations (47) through (50) describe 
both convectively stable and convectively unstable regions, but D and i/diss are effectively 
zero in stable regions that are not very close to marginal stability. 

A few words on the relationship between the two-fluid convection model and standard 
mixing length theory are in order. The left-hand side of equation (47) is reminiscent of 
the divergence of the convective heat flux that appears in standard mixing length theory 
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and discussions of forced turbulent mixing in clusters. However, the convective internal- 
energy flux (SvrSp)/^ — 1) on the left-hand side of equation (47) is not proportional to the 
specific-entropy gradient. In fact, (Sp Sv r ) vanishes as p CT — > 0, because the total pressure 
perturbation in a fluid element is assumed to vanish. However, the approach taken in this 
paper would be similar to standard mixing length theory in the p cr — > limit if the average 
velocity were retained. In a convectively unstable region, downward falling fluid elements 
are denser on average than upwardly rising fluid elements. The random velocities are thus 
associated with an inward mass flux. In standard mixing length theory, the total mass flux 
vanishes because a steady state without mass sources or sinks is assumed (Cox & Giugli 
1968). This requires a non- vanishing outward average radial velocity, 

M = - (51) 

If we add equation (36) to the dot product of v with the steady-state momentum equation, 
employ V • (pv) = 0, and take the limit p cr — > 0, we obtain the total-energy equation 

= -V- (^ + ^ + P v$ + *v) +H tc -R, (52) 

where 7pv/(7 — 1) is the enthalpy flux, and F v is the viscous energy flux [see, e.g., Rudiger 
(1989), equation 8.40]. If one sets 

(Sp Sv r ) = c mix Ap sgn(Z)w, (53) 

the enthalpy flux associated with the average radial velocity of equation (51) is 

iMM = _ DpT d ° (54) 
7 — 1 dr 

where the angle brackets around average quantities have been dropped on the right-hand 
side, s = Cy ln(p/p 7 ) is the specific entropy, and Cy = e/T is the specific heat at constant 
volume. The enthalpy flux in equation (54) is the same as the heat flux in standard mixing 
length theory. 

Although (v r ) must be included to recover standard mixing length theory in the p cr — > 
limit, (v r ) can be safely neglected in the energy equation when the model is applied to 
clusters. In the limit of vanishing M, the (v r )(p) terms in the average of equation (36) are 
smaller than the (5v r 5p) terms, which can be seen as follows. Using equations (51) and (53), 
one can write 

(v r )(p) _ Ap (Apy 1 

(5v r 5p) ~ p\p) ■ (55) 
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This ratio is small for the following reasons. Convective regions must be close to marginal 
convective stability to produce the subsonic turbulent velocities required to balance cooling, 
which implies |Ap| <C \lp/r\. This condition, coupled with equations (25) and (26), gives 

^ (56) 
Cy or 

Equation (56) implies that |Ap| ~ Ip/r, since ds/dr is observed to be ~ Cy/r in clusters. 
Thus, |Ap/p| ■ \Apfp\~ 1 <C l. 4 It is true that if M is sufficiently large, then the (v r ) terms 
in the energy equation are comparable to the radiative-cooling term, as in the cooling flow 
model. However, in the model solutions to be presented below, M is a factor 10 2 — 10 3 
smaller than in the cooling-flow model, so that the neglect of (v r ) is okay. 5 



Equation (42) can be rewritten 



Q = -Ptot-r- + 
dr 



(7 ~ 7cr)XPc 
IP + IcrPcr 



dp 



tot 



dr 



(57) 



where 



X = — ■ (58) 

Ptot 

One expects dx/dr > throughout most of a cluster since cosmic rays contribute a larger 
fraction of the pressure nearer the central radio source. Since 7 > 7 cr , one expects Q to be 
negative in clusters. Thus, if a fluid parcel is displaced inwards (/ < 0), the thermal pressure 
inside the displaced parcel is larger than in the surrounding medium (IQ > 0), essentially 
because the total pressure is the same inside and outside the parcel and x is larger inside. 
Since the inwardly displaced parcel gets compressed, the first term on the right-hand side of 
equation (47) (— (opV • ov)) is positive, and acts as an energy source for the thermal plasma. 
This source term can be interpreted as in part due to pdV work on the thermal plasma by 
the cosmic-ray pressure, which helps compress the displaced parcel. 

Equations (47) through (49) form a system of two second-order equations and one first- 
order equation for p, T, and p cr . Five boundary conditions are required to specify a solution. 
Two boundary conditions are obtained by imposing a density p uter and temperature T outer at 



4 If one takes pg ~ p/r and V • (<5v Sp) ~ (Sv r 5p)/r, then f/diss/V • (5v dp) ~ (Ap /p)(Ap /p)^ 1 , which is 
again small because convective regions are near marginal convective stability. Thus, dissipation of turbulent 
motions is less important than turbulent mixing in the model of this paper, as is borne out by numerical 
solutions to equations (47) through (49). On the other hand, if an external source of turbulent motions is 
invoked as in, e.g., El-Zant et al 2004, then viscous dissipation may be more important in comparison to 
turbulent mixing than in the present model. 

5 Very close to cluster center, (v r ) becomes large even for small M, and the model is inaccurate. 
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radius r outer . Two additional boundary conditions are obtained by taking dT/dr and dp cr /dr 
to vanish at the origin. The fifth condition is obtained by assuming that p CT — > as r — > oo. 
This condition is translated into an approximate condition on p CT at r outer as follows. The 
value of r outcr is chosen to be much greater than r s and D Q /vd, so that for r > r outcr , S* ~ and 
D cr ~ t^r. Equation (48) then implies that p CT ~ C\ + C2r~ 2 for r > r outer assuming the ICM 
is convectively stable and far from marginal stability. Since p CT — > as r — > oo, ci = 0. The 
fifth boundary condition is then taken to be dp CT /dr = —2p cr /r at r outer . Numerical solutions 
are obtained using a shooting method. Values are guessed for p, T, and p cr at r = and 
the equations are integrated from to r out er- The guesses are then updated using Newton's 
method until the three boundary conditions at r outcr are met. 

Figure 1 shows an example solution for the following parameters: M bh = 1O 9 M , rj = 
0.003, 7 = 5/3, 7 cr = 4/3, a — 1, u' = 0.1 km s _1 kpc _1 , r s = 25 kpc, v d = 18 km/s, 
-D = 10 28 cm 2 /s, and 9 = 1/3. The values Cdi SS = 0.42 and c mix = 0.11 are adopted based on 
a number of previous studies, as discussed by Dennis & Chandran (2004). The gravitational 
potential is taken to be 
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with v c = 1150 km/s, r c = 40 kpc, v g = 400 km/s, and r g = 1 kpc. The term proportional 
to v 2 (v 2 ) represents the cluster (central-galaxy) potential. A mean molecular weight per 
electron p/n e m H of 1.18 is assumed, where m# is the mass of a hydrogen atom, and the 
ratio ni/n e is set equal to 0.9. 6 The values p ou ter = 5.38 x 10~ 3 cm" 3 and T outcr = 7.91 keV 
at r outcr = 200 kpc are determined by linear interpolation between data points for Abell 478. 
Density and temperature data for Abell 478 from Chandra observations, provided by S. Allen, 
are also plotted in figure 1. Abell 478 is a highly relaxed, x-ray luminous cluster, with a 
cooling-flow-model mass accretion rate of ~ 1OOOM yr" 1 , and a comparatively weak central 
radio source (Allen et al 1993, Sun et al 2003). 

It can be seen from the figure that the the model density is too small for r < 50 kpc. This 
discrepancy arises because the specific entropy of the thermal plasma increases outwards, and 
thus the cosmic-ray pressure gradient must be large in order for the ICM to be convectively 
unstable and for turbulent heating to help balance radiative cooling. When cosmic rays 
contribute a large fraction of the total pressure support, the thermal pressure gradient does 
not have to be as large, and thus the central plasma density does not have to be as large. A 
variety of model parameters have been investigated with similar results. 



6 These values are appropriate for a hydrogen fraction X = 0.7 and a helium fraction Y = 0.28, with fully 
ionized hydrogen and helium (Zakamska & Narayan 2003). 
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Fig. 1. — Top panels: comparison of model density and temperature with observations of Abell 478. Data 
points and one-sigma error bars (provided by S. Allen — see Sun et al 2003) assume Qm = 0.3, Oa = 0.7, 
and Hq = 70 km s _1 Mpc _1 . Bottom panels give the ratio of cosmic-ray pressure to thermal pressure and 
the turbulent velocity as a function of distance from cluster center. 

The turbulent velocity in figure 1 is also surprisingly small. This is in part because 
the model underestimates the plasma density in the central 50 kpc, which weakens radiative 
cooling and the need for turbulent heating. Figure 2 shows a numerical solution for the 
parameters v c = 1500 km/s, x c = 10 kpc, v g — 0, 9 — 0.37, Vd = 20 km/s, and D = 
3 x 10 28 cm 2 /s, with other parameters the same as in figure 1. The enhanced gravitational 
acceleration increases the density and radiative cooling, which results in the need for greater 
turbulent heating and a larger u(r). These parameters, however, are unrealistic since they 
imply a velocity dispersion of ~ 1000 km/s within the central galaxy at r = 10 kpc, which 
is much larger than observed velocity dispersions in giant elliptical galaxies. 

A further difficulty with the model is that a critical point is encountered for sufficiently 
large Vd and/or D , or for sufficiently small 9. It can be shown that a sufficient condition 
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kpc kpc 

Fig. 2. — Same as figure 1, but with different model parameters, as described in the text. 

for avoiding a critical point is that (7 cr — l)D cr /(7 — 1) is everywhere less than K T /n e kB, 
the diffusion coefficient of heat-carrying electrons, and the models plotted in figures 1 and 2 
satisfy this inequality. The potential for a critical point raises a number of interesting 
questions, but these lie beyond the scope of this paper. 

For reference, in the numerical solution plotted in figure 1, the Bondi accretion rate is 
1.2M yr" 1 , L cr = 2.1 x 10 44 ergs/s, the radiative luminosity out to 200 kpc is 1.9 x 10 45 erg/s, 
and the radiative luminosity out to 600 kpc is 4.9 x 10 45 ergs/s. For the model plotted in 
figure 2, the Bondi accretion rate is 5.0M Q yr -1 , L CT = 8.5 x 10 44 erg/s, the radiative 
luminosity out to 200 kpc is 2.2 x 10 45 erg/s, and the radiative luminosity out to 600 kpc is 
4.5 x 10 45 erg/s. (The luminosity out to 600 kpc is smaller in the model of figure 2 because 
the larger gravitational acceleration causes the density to drop off faster at large r.) 
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3. Summary 

AGN feedback is a promising explanation for the heating of intracluster plasmas on 
both observational and theoretical grounds: there is an active radio source at the center of 
almost every strongly cooling cluster (Eilek 2003), and AGN feedback in combination with 
thermal conduction can lead to globally stable equilibria (Rosner & Tucker 1989, Ruszkowski 
& Begelman 2002). The recent detection of moderately subsonic plasma motions in the 
Perseus cluster (Churazov et al 2004) further suggests that turbulent mixing and/or turbulent 
dissipation plays an important role in heating the intracluster medium. This observation also 
poses a theoretical challenge. For velocities in the range of those observed, turbulent heating 
can either overwhelm radiative cooling or be too small to offset cooling, depending on the 
precise value of the rms velocity u and the velocity length scale. If turbulent heating indeed 
balances radiative cooling, an explanation is needed for how u is fine-tuned to the required 
value. 

The model presented in this paper connects AGN feedback to turbulent heating, and 
seeks to explain the fine-tuning of u. In the model, radiative cooling of intracluster plasma 
is balanced by a combination of turbulent heating and thermal conduction. The turbulence 
is generated by convection, which in turn is produced by the buoyancy of cosmic rays gen- 
erated by a central radio source. Convection is a natural way for a central relativistic jet to 
generate turbulence due to the jet's small momentum-flux-to-energy-flux ratio: rather than a 
relativistic outflow stirring intracluster plasma like an oar in water, the radio source inflates 
the central region, allowing gravity to generate plasma momentum after parcels of plasma 
and cosmic rays in the central region become buoyant. A two-fluid (plasma and cosmic ray) 
mixing length theory is developed to treat convection in the ICM. By linking the turbulence 
amplitude to the mass accretion rate and AGN feedback luminosity, the model provides an 
explanation for how the turbulent velocity achieves the value required for heating to bal- 
ance cooling. Although stability is not investigated in this paper, it is expected that the 
model equilibria are stable for the same reasons as other AGN- feedback/thermal-conduction 
models. 

The model provides equilibrium solutions for n e (r), T(r), p cr (r), and u(r) once the 
plasma temperature and density are specified at some suitably large outer radius. For typical 
cluster parameters, the ICM stays fairly close to marginal convective stability in convective 
regions, implying subsonic turbulent velocities. The principal shortcoming of the model is 
that the model density is too small within the central ~ 50 kpc, where the model ICM is 
convective. The reason is that the cosmic-ray pressure gradient must be large for cosmic- 
ray buoyancy to generate convection. This in turn reduces the thermal pressure gradient 
in the central region and decreases the central plasma density. A more realistic model of 
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convection in the ICM in which relativistic and thermal plasmas are only partially mixed 
may have greater success. 

I thank Steve Allen for providing the data for Abell 478, Eric Blackman, Eliot Quataert, 
and the anonymous referee for providing helpful comments on an earlier version of this paper, 
and Steve Allen, Steve Cowley, Ue-li Pen, and Steve Spangler for helpful discussions. This 
work was supported by NSF grant AST-0098086 and DOE grants DE-FG02-01ER54658 and 
DE-FC02-01ER54651 at the University of Iowa. 
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